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Abstract 

In ten dimensional type II superstring, all perturbative massive states are unstable, 
typically with a short lifetime compared to the string scale. We find that the lifetime 
of the average string state of mass M has the asymptotic form T < const. g~ 2 A/ -1 . 
The most stable string state seems to be a certain state with high angular momentum 
which can be classically viewed as a circular string rotating in several planes ( "the rotating 
ring"), predominantly decaying by radiating soft massless NS-NS particles, with a lifetime 
T = cog~ 2 M 5 . Remarkably, the dominant channel is the decay into a similar rotating ring 
state of smaller mass. The total lifetime to shrink to zero size is ~ M . In the presence 
of D branes, decay channels involving open strings in the final state are exponentially 
suppressed, so the lifetime is still proportional to M 5 , except for a D brane at a special 
angle or flux. For large mass, the spectrum for massless emission exhibits qualitative 
features typical of a thermal spectrum, such as a maximum and an exponential tail. We 
also discuss the decay properties of rotating rings in the case of compact dimensions. 
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1 Introduction 

Understanding the decay properties of massive string states may uncover new aspects of string 
theory dynamics, and it may also have direct implications in cosmology (see e.g. [1-4]). Ac- 
cording to their decay properties, the quantum states of the type II string spectrum can be 
divided into two classes: 

i) Those for which the dominant decay channel is by massless emission. 

ii) Those for which the dominant decay channel is by emission of massive particles. 

In [5] we argued that massless emission is the dominant channel whenever the quantum 
state corresponds to large closed strings which cannot break during the evolution. An example 
of this is the rotating ring: the string is circular, and there is never a contact between two 
points of the string. For these class i) string states, the decay channels into massive states are 
exponentially suppressed, we indeed find r massivc = 0(e~ cM ), where the value of c depends on 
the specific state. 1 An example of a class ii) state is the folded closed string corresponding to 
the quantum state with maximum angular momentum: the points of the string are always in 
contact with the other side so the string can break at any moment. This picture was confirmed 
in detail by the calculations of [5]. With this line of argument, in [5] we have identified a 
one-parameter family of long-lived string states. 



1 The numerical constant c > can be very low if there are two points of the string that are close to each 
other. In such a case the decay into massive states can be viewed scmiclassically as tunnelling effect. 
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For the class i) quantum states, the lifetime is given by 7~W = r~^, sless , while for the class 
ii) quantum states, the decay rate for massless emission r masslcss puts an upper bound on the 
lifetime, i.e. < r~^ ssless . In this paper we will concentrate on the calculation of r massless , 
first finding a simple general formula, and then carrying out the analytic calculation in detail 
for the case of the rotating ring. The analytic calculation confirms the lifetime % mg = p-M 5 , in 
ten uncompact spacetime dimensions, obtained in [5] by numerical evaluation of the imaginary 
part of the one- loop correction to the two-point function (here we also determine cq). Further, 
we show that for large M the decay occurs by the emission of massless quanta of feeble energy 
oj ~ 1/M and that the ring final state is just the same as the initial one, only slightly shrunk to 
conserve energy. Moreover, we show that for large mass a classical computation gives exactly 
the same results as the quantum one. That is, the rotating ring slowly radiates as a classical 
antenna (section 3.3). In section 4 we show that the rotating ring has the same lifetime 
%mg = -%M 5 in the presence of generic configurations of D branes. This is due to the fact that 
the decay into open strings is exponentially suppressed, except for a D brane at a special angle. 

The case of compact dimensions is discussed in section 5. The decay properties of the ring 
then depend on the ring configuration, in particular, on whether the ring lies on compact or 
uncompact dimensions. For configurations where decay into winding modes is suppressed, the 
lifetime of the rotating ring is T = const. g~ 2 M d ~ 4 , where d is the number of uncompact 
spatial dimensions, while the total lifetime for the ring to shrink to microscopic size is T tot = 
const. g~ 2 M d ~ 2 . 

In section 6, we compute f mass i ess representing the average over all initial states of given 
mass M (following the similar calculation carried out in [9] for the bosonic string). 

Some earlier discussions of decays of massive string states can be found in [6-8]. The 
average decay rate by massless emission in bosonic string theory is given in [9] and a study of 
massive emission from an average string state is in [10] (averages involving two massless vertices 
were discussed in [11]). The classical breaking of the string and the correspondence with the 
quantum calculation is discussed in [12]. A numerical study of the decay of the rotating ring 
in compactified dimensions from the one-loop two-point function is in [13]. 

2 General formula for decay rates 

The total decay amplitude of a state |$o) of mass M for emission of a massless (fermion or 
boson) particle (oj, k) represented by a vertex operator V is given by 

= r dn^duj u d - 2 _ E ,_ u) = 9\ ^ d - 2 r 

y J 8M J (2n) d - l E' v 8(27r) d - 1 M 2 J H ' v ; 

where d is the number of spatial dimensions and 

a R = J2(®o\V R (k, iy\$ N ,)($ N .\V R (k, 1)|$ ) , (2.2) 

*JV' 

and a similar expression for ul- The sum goes over all states with mass M' = vW 7 , with energy 
E' = a/M /2 + oj 2 , and (in units a' = 4) 

,= M2 - M/2 =^L, M 2 = N, N = N-N>. 
2M 2VN 
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Introducing a projector onto states of level N', 



£|<M<<M = ^>- 7V ', (2.3) 



the decay amplitude can be computed as 



a R =<b- z- N '(%\V R (k, lYz^V R (k, 1)|$ ) . (2.4) 



z 

By using the standard commutators 

[a, a f ] = 1 , V 1 } = 1 , ^ nata e a = e fl/2 "z mtfl , z^e* = e^ /zT z ri>H , 
one can write 

z**V R (k,l)=V R (k,z)z** . (2.5) 

Thus 

a R = z N - N '(<f>o\V R (k,iyV R (k,z)\<!> ) . (2.6) 

For massless NS-NS emission, one considers the vertex operator V = g s VlV r , with 

V R = (£.dX R + % ^ R k.jp R ) e tk - XR , V L = (ldX L + % i.jp L k.^ L ) e lk - x \ (2.7) 

f .Jfe = f.Jfe = , k 2 = . (2.8) 

We will work in the gauge where £ = £o = 0. Equations 1)2. ip . ()2.6|) . allow to compute 
the decay rate of any quantum state of the string spectrum by computing a matrix element 
involving standard free/creation annihilation operators. In the next section we carry out this 
calculation for the rotating ring. The calculation can also be carried out for other states in a 
straightforward way. 



3 The rotating ring 

3.1 Decay rate 

We now consider the following state 

where we define = b\ + = a ~ 1 J^' 1 , $ = c\ + = J^^ 1 (see appendix A for notation). The 
operator b> (£•) adds one unit of Right (Left) angular momentum in the plane 12 (34). This 
state has 

N R = N L = N , a'M 2 = AN , (3.2) 
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Jl2R = N, Ji 2L = 0, J 3 4R = 0, J 3 4L = iV, (3.3) 

so that Ji 2 = J 34 = N. The classical string soliton solution with the same values of J\ 2 r, Jyil, 
J34R, J34L is as follows 



Zi — L e"*^) , Z 2 = L , L=^^, (3.4) 

X = 2v^L r , a e [0, 2tt) . 

It describes a circular string rotating around its axes in two orthogonal planes. One can also 
write similar rotating ring solutions (and the corresponding quantum states) with rotation in 
different planes. 

The lifetime of this state was studied in [5] by a numerical evaluation of the imaginary part 
of the one- loop two point amplitude. It was shown that it goes like T ~ M 5 , and that the 
dominant decay channel is massless emission. Other decay channels representing decay into 
two massive particles are suppressed exponentially for large mass of the initial state, i.e. the 
decay rate for these massive channels is T ~ exp(— cM 2 ). Thus in the large mass limit, the only 
contribution which is not exponentially suppressed is decay by emission of a single massless 
particle (multiple massless particle emission is suppressed by extra factors of g 2 ). This can be 
computed analytically by using the operator approach, as explained in section 2. 

Let us start considering emission of a NS-NS massless particle. It is convenient to introduce 
complex notation 

k = —j={k\ + ik 2 ) , k=—(k 3 + ik 4 ), (3.5) 
v2 v 2 

£ =-L(£ 1 + i£ 2 ) , i= -L (f 3 + zf 4 ) . (3.6) 

Now consider the computation of <jr ()2.6|) . 

r Hz h N - tf N 

°* = f™ ^H^^^/2V R (k^V R (k,z)^^ 1/2 \0) , (3.7) 

with N = N — N'. The product Vn(iyVji(z) gives rise to five contributions, 

<?r = v\ + o"2 + 0-3 + 04 + 05 , (3.8) 

with 

ai = __L j *E z No (0\b N : e- lk - xil) e ik - x{z) i.dX{l)i.dX{z) : tf N \0) , (3.9) 



where we have used (O^f tp_ 1 10) = 1, and <t 2 , ■■■,<j 5 are given in appendix B. 

2 2 

For large N, <j\ turns out to be the dominant contribution, while <t 2 , ...,cr 5 are suppressed 
by inverse powers of N (see Section 3.2). 

Using the algebra of creation/annihilation operators, we obtain o\ = a[ r ^ + cr[ H \ with 



a, -p+^jf / ^m\{m + 2)\{N-m-2)\ Z ™ 



m=0 



in) _ oee I dz z N {z + £ iV!(2^)- (l-zf 

m=0 



m! 2 (iV — m — 1)! z r 



Now we expand the binomial (1 — z) r and perform the integrals in z around a small contour 
near z = 0. We find 

N-No-l 

a{ n = A(2 K R) N ^ 1 (e 2 n 2 + e 2 R 2 ) £ C/ (r) (-2^) r , (3.10) 

r=0 



JV-jV 

<jj 77) = A(2^R) N °~ l ee ^ c 77 (r) (-2^) r , (3.11) 

r=0 

where 

C/W = jV!(2r + 2JV )! 

V; r\(N - N Q - 1 - r)l(r + 2N )l(r + N - l)\(r + N Q + 1)\ ' V ' 7 

iV!(2r + 2iV - 2)![r 2 + (r + iV )(2iV - 1)] 
C/Arj r!(iV-Ar -r)!(r + 2iVo)!(r + iVo-l)! 2 ' 1 1 

Next, we sum over polarizations, using the relations given in appendix A. We obtain 

ffi> = - 2 (^) S c,(r) \^r) • (3 ' 14) 



where v 2 = sin 2 # cos 2 a cos 2 /5. The Left contributions are similar, changing k, Rbj R,R and £ M 
by I'*, with 2k 2 = J§ e 2# sin 2 9 cos 2 a sin 2 (3. 

Now we consider the product <jt = ®l x cr and perform the integral over the angular part. 
Defining 

J(r, s ) = [ d 8 tt(sm 6 cos a) 2r+2s (cos f3) 2r (sin f3) 2s = ( 2?r ) 4 v 7 ^ r\s\ (g lg) 

y 8 r(| + r + s) 



we find 



y rf 8 fi a T = E n + £ 22 + 2£ 12 , (3.17) 



with 

r2 \ 2V JV-JVo-1 / at2 \ r+s 



^ 11=4 \4n) E [mj ai^iWir + No + hs + No + l) , (3.18) 



r,s=0 
iV-iVo 

= 16 (iA7) E (-4^) cjr(r) C/J r( a )J(r + JVb-l,a + JVb-l), (3.19) 

/ r,s=0 ^ ' 



MY 



/ at2 \ 2N ~1 N-N -l N-N , N 2 \ r+s 

2 = "8 ^ £ E ( "4^ ) c 7 (r)c„( a )iif (r + 7V - 1, S + 7V - 1) , (3.20) 



111 

J(m, n) = I(m, n) — -I(m + 1, n) — -I(m, n + 1) + ~^I\ m + l,n + 1) , 

K(m, n) = J(m + 2, n) — -I(m + 2, n + 1) . 

Thus the total decay amplitude for emission of a massless NS-NS particle from the rotating 
ring is given by 

rM = ^^(s 11 + s 22 + 2S 12 ) . (3.21) 



Using u = A r o/(2M), m + v AT, this can be written as 

r(u)=g 2 s N-l f{N ,N) , (3.22) 

with 

n 2 N 7 

/(^,^) = ^ ? ^|(Sh + S 22 + 2E 12 ) . (3.23) 

This function f(N , N) becomes independent of iV for large N. This is due to the fact that 
the main contributions arise for finite N for any A" (see sect. 3.4). In fact, it has a maximum 
at iVo = 3 irrespective of the value of iV (provided TV > 9). For large N, one can then use the 
Stirling approximation to write, in cj(r), cu(r), 

N\ 

— ~ N m . (3.24) 

(N-m)\ 

Then the A" dependence in f(No, N) fully cancels out. Indeed, from ci(r)ci(s) one gets a power 
jy2N +r+s+2 cance ls all powers of A" in En, and similarly for Ei 2 and S 22 . Thus we can 

replace f(N , N) by f(N ) = f(N , 00) so that the large A" result is 

r(AT ) = g 2 M- 5 f(N ) . (3.25) 

The total decay rate can now be computed by summing over all N from 1 to A^: 

N 00 

r mass iess = g 2 s M" 5 Wo, N) = g 2 s M" 5 ^ f(N ) , (3.26) 

N =l N =l 

In the next section, we show that massless fermion emission and emission of massless RR states 
are suppressed by factors 1/N and 1/N 2 respectively, so r mass i ess represents the complete decay 
rate modulo 1/N corrections. Thus, for large initial mass, the lifetime of the rotating ring is 
given by 

^ ing = r-L lcss = c s; 2 M 5 , c ee f; f(N ) - . (3.27) 

N =l 

This exactly reproduces the result of [5] obtained by computing the imaginary part of the one- 
loop two-point function (modulo the numerical coefficient Co which had not been determined 
in 
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3.2 Dominant decay channel 

The large TV behavior of the various terms giving the massless (radiation) emission rate can be 
understood in brief by the following sample computation. 

Consider the term ($™ g |0{X}|$™ lg ), 6{X} =: ^.dX{l)tdX{z)e- ikX ^ +ikX ^ : , appear- 
ing in the Right factor. We can restrict the full operator 0{X} to the part containing b, b\ as 
only this part gives a contribution. Ignoring the factors £ and the dependence on z, we have 

<•&■ \6{X}\^) ~ <0|^=6te-^'6 e«* ^=|0> = £ . 

* V • m=0 

For large AT, one can use Stirling formula (jH.24)) . One obtains 

(^ s |OWI$7> - N J2 { ~ 2KR * )m = Nfi(N ) . (3.28) 

We have used 2kR = k\ + fcf = u; 2 f 2 = Thus this matrix element scales like iV for large 

N. 

Another term (called a§ in appendix B) comes from the contraction 

(CdX(l)£.dX(z))($i ng \ : e -ik-X(l)+ik-X(z) . | $ ring^ _ ^ ^ 

This gives rise to the following factor: 

<°'S 5' 0> = t §^ - fM) ■ (3.30) 

v v • m=0 

which is of order 1 at large N. Therefore it is subleading by a factor 1/N with respect to 
expression ()3.28j) . The origin of the extra factor of N in (I3.28|) can be traced back to the action 
of dXdX on the ring state \&n)- in (13.29)1 . the operators dXdX are contracted and they do 
not give A" dependence. 

Similarly, the contribution from the term in the vertex (£ • -0 k k -ip) (computed in 

appendix B) will be subleading for a factor 1/N 2 , since there is no dXdX operator and there 
is an extra power 1/N coming from ~ l/\f~N. 

For the same reason, decays by massless fermion and RR emission are suppressed. Indeed, 
in the picture we are working, the vertex operator for a massless fermion in the RNS sector is 
given by V f = V£V L , with 

Vl{u a ; fc") = e&- u a Q a e ikXR , (3.31) 
& a = D 1 , 1 ...Dl 1 , -D?i =e ± ^ m 

-*-2 -"-2 - c 2 



where tp and m represent, respectively, the bosonized /?7 system and the bosonized ip^. Here 
Vl is the NS Left vertex operator (see (|2.7|) ). 

In view of the above discussion, since for these vertices there is no factor £,.dX, emission 
of V* is suppressed by a power of order 1/N with respect to emission of the NSNS massless 
particle (j2.7|) . and emission of a massless RR state V* = V^Vf is suppressed by a power 1/N 2 . 
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We have checked this suppression by a numerical evaluation of massless fermion and RR 
emission. The calculation follows [5] , based on evaluating the emission rate from the imaginary 
part of the correlator of two vertices on the torus, and looking at the contributions from the 
different spin structures. We have found that RNS and RR are subleading, by factors 1/N and 
1/N 2 respectively, with respect to NSNS, as expected in view of the above argument. 

The previous considerations lead to an important consequence: for large initial mass the 
dominant decay channel is that where the final state is a rotating ring of lower mass, i.e. similar 
to the initial state with the replacement iV — > N — Nq. The leading process is 

®7 g - + 7H , (3-32) 

where 7 is the massless particle. To see this explicitly, we start with the dominant contribution 
(|3.9jl . As pointed out above, because of the normal ordering, we can drop all other mode 

operators in X, keeping only b and b\ Defining Y(z) = X(z) we write 

6,6t 



~<f— z No (0\b N : e- ikY ^i.dY{l)e ikY ^i.dY{z) : tf N \0) . (3.33) 



Next, using the above result (see also appendix B) that the contraction of dX(l)dX(z) gives 
rise to a contribution suppressed by a power 1/N, we can write 



. e ~ik.Y(l) f 



^.dY(l)e ikY{z) £.dY(z) :=: e~ ikY(1) CdY(l) : : e ikY{z) £.dY(z) : M +0(-^ ■ 
Therefore o\ can also be written as follows: 

ai = -Jfl X>l 6 " ■ e- ikY{1 \.dY(l) : \* N .)(* N ,\ : e tkY ^^.dY(l) : 6^|0) . 

So consider the calculation of ($n'\ '■ e lkY ^C,.dY(l) : tf N \0) . The operators e lkYiY) and dY(l) 
contain only b, 6' oscillators, Therefore, the states |$jv'} can only be made by applying 
to the vacuum. But there is a single state with Nr\$n') = N'\$n') made with b\ namely 
l^iv) = (b ) N |0). Thus in the large limit the leading contribution to the decay rate in the 
sum over all \<&n') comes from a circular ring of lower mass M' = yN'. 

As an independent check of the decay rate of section 3.1, we now compute the matrix 
element 

b N> tf N 
A R = (0|-= : {.dY(l)e ikY M : | ) , (3.34) 



'N'\ 

corresponding to the leading process (J3.32j) . In the large N limit, this should coincide with the 
decay rate of section 3.1 obtained by summing over all final states of mass M' = y/N'. We find 



A 



s=0 



s\(N- N -s)\ 



en(N-N -s) e 1 
_ (A^o + s + 1)! + 2^(N Q + s- 1)!_ 



(3.35) 
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For large N, using the Stirling formula (|3.24j) this becomes 



A R = -2iVN (-V2Nk) No J2 



-2kkNY 



s=Q 



Nek 



(N + s + l)\ 2x/Nk(N + s-1)\ 



(3.36) 



Since k is proportional to 1/yN, the full expression is proportional to y/N. Now consider 
|t4#| 2 . This can be written as 

\A R \ 2 = aW + aM , (3.37) 



with 



=AN 2 (2KRN) Na -\e 2 K 2 + e 2 R 2 ) V -— ( 2kRN J +S ( 3 .38) 

' ^ r\s\(N + s + l)!(iVo + r - 1)! 



r,s=0 



= 2Nee(2kRN\ N °- 1 V (~ 2 ^) r+ K 2 ^) 2 + (No + s)(N + s+l)(N +r)(N + r + l)] 
1 ' r!s!(iV + s + l)!(iVo + r + l)! 

(3.39) 

By defining a new summation variable r' = r + s, and performing the sum over r at fixed r', we 
find that exactly reproduce a[ J \ a[ H ^ given in eqs. (|3.10jh (|3.11|) . with cj(r), cu(r) 

given by the Stirling limit ()3.24|) of (|3.12|) , (|3.13|) . This provides another derivation of the decay 
rate of the rotating ring T = const M 5 , by computing only the dominant decay channel. 

3.3 Rotating ring as a classical antenna 

Having shown that the dominant decay channel of the rotating ring is by massless emission, it 
is instructive to see that the exact large mass decay rate can be reproduced by regarding the 
rotating ring as a classical antenna. 

The classical radiation from a source T^ u (x ,x) is given by [14] 

a 

J = J dx d d x e -*°-^ T^xq, x) . 
For a classical string solution, the energy momentum tensor is 

Tp, = J dadr 5 {d) (x - X(t, a))d^d+X v . (3.41) 

For the rotating ring ()3.4j) . 

J = J dadTe 2iNoT - it ^ T) ^td-X lx d + X v , (3.42) 

ik.X(a, t) = ikZ\ + %kZ\ + iRZ 2 + ikZ\ , 
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where Z l5 Z 2 are given in (|3.4j) . This can be written as J = Jr Jl, with 

J R = da_ e lNaa -- mLe ^~ iKLe ^~ (ieLe ia - -ieLe~ lr7 -) , (3.43) 
Jo 

r-2-K 

J L = da + e iNoa+ - fkLeia+ - ikLe ^ + {iiLe i(7+ -ieLe- icT+ ) , (3.44) 
Jo 

with o± = t ± a and L = y/2N. Performing the integrals over a + and cr_,we obtain 

Jr = 2^-^)- 1 j(-2^-( (Ab ;; + 1} , + 2JVK(jV ; + , _ ■ <^> 

Remarkably, this coincides exactly with ()3.36|) (up to an irrelevant phase i N ° which drops out in 
| J| 2 ). It is important to note is that the quantum calculation reduces to the classical radiation 
process after taking the Stirling limit (J3.24j) . i.e. for large mass M or large excitation number 
N. This is the limit where a classical description in terms of a string soliton can be applied. 



3.4 Spectrum of massless emission 

The spectrum of massless emission M 5 x T(u) is shown in fig. 1 for different values N = 
20, 50, 100. A salient feature of this spectrum is that the maximum is always at iV"o = N—N' = 3 
for all iV > 9. The large N form of the spectrum is given by (|3.25|) . which in terms of uo reads 

T(u) = g 2 s M- 5 f(2Mu) . 

Thus the rotating ring decays mainly by emitting soft NS-NS massless particles (gravitons and 
Bp,) of frequencies tu = N /(2M) ~ M" 1 . 



1. -10 



5. -10 



I l 



10 12 



Figure 1: Spectrum M 5 x T(u) of massless emission from the rotating ring as a function of 
iV = 2Mu, for N = M 2 = 20 (lower points), N = 50 and N = 100 (higher points). Here 
g 2 s = 0.1 and a' = 4. 

Since the dominant channel is the one where the final state is a rotating ring of lower mass, 
the decay can be pictured as a rotating circular string whose radius gradually decreases as it 
loses energy by radiation. 
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For large N Q , the spectrum exponentially goes to zero. It is of the form (see appendix C) 

r(w) = N~l e -2JVo(21og2-l) = N -l e -% ? (3 _ 46) 

To = S' «o 1 = 4(21og2-l) . 

The exponential form of the spectrum at large u, and the presence of a maximum at finite 
frequency allows one to approximate the spectrum curve by a thermal spectrum by introducing 
a "grey body factor" cr grcy (c<j) defined by 

r(u) = const. M~ 5 u d ~ 1 a grcy (uj, T ) — 3- , 

e T o - 1 

i.e. 

a gvcy (u,T ) = M 5 r(^)^ 1 - d (e^ - 1) . 

where T(u) is the expression obtained in subsection 3.1. The grey body factor here is introduced 
of course artificially, but given the decay properties of the rotating ring, as a smooth radiation 
process with frequencies peaked at u ~ M _1 , it is useful to picture it as a thermal radiation 
process, even if the exact spectrum formula is actually more complicated. 

We have seen that the spectrum can be also obtained by a classical computation; the fact 
that it is peaked at small frequencies is in agreement with the classical analysis of [3] in the 
case of a smooth antenna. 

We have also examined the spectrum of the string state with maximum angular momentum 
(6'6') iV |0). In this case we find that the spectrum does not resemble a thermal spectrum, in 
particular, it is not exponentially suppressed at large frequencies (more precisely, in the region 
1< iV < N). 



4 Presence of D branes 

A question of interest is what is the lifetime of the rotating ring in the presence of D branes. 
Here we show that generically the ring cannot break into open string fragments, with Neumann 
or Dirichlet boundary conditions, nor in the presence of generic fluxes. 

Moreover the rotating ring cannot either decay into another closed string plus open string 
fragments. All these processes are exponentially suppressed, i.e they are of order 0(e~ M ). 

There is an exceptional case in which the breaking is not exponentially suppressed, namely 
if there is a D brane tilted with an angle equal to ±7r/4 with respect to the Z\ and the Z 2 
planes. This configuration is equivalent to a D brane such that its projection in the Z\,Z 2 
subspace appears as a D4 brane with a suitable flux. 

We begin by discussing the picture by the classical world-sheet equations. The classical so- 
lution corresponding to the circular string is given in eq. (|3.4j) . which we repeat for convenience: 

X^Gjir-a), X c k = G k {r + a), (4.1) 

where j — 1,2 and k = 3,4, with G 1:3 (x) = L cos(x), £2,4 0*0 — L sin(x). 

The classical equations are second order in r therefore the breaking must require continuity 
in r of the classical functions and of their first derivative (see [12]). In particular, this initial 
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condition ensures that the Virasoro constraints are satisfied for the outgoing fragments. This 
initial condition holds for every coordinate. 

Take generically an open string to be of the form X open (r, a) = f(r + a) + g(r — a). For 
M Neumann or D Dirichlet boundary condition g(x) = ±/(x), but we consider more general 
cases resulting from rotations or presence of fluxes. 

Suppose that at time r = the closed string breaks and it becomes an open string. The 
initial condition requires 

Z c (0,a) = Z^ n (0,a) = f(a)+g(-a), (4.2) 
d T Z c (0, a) = Td„Z c (0, a) = d T Z open (0, a) = f(a) + g'(-a) . (4.3) 

with =F = — for the j and = + for the k components. Since the continuity requirement (|4.2jl 
must hold for every a, we can take the derivative with respect to a and get 

d a Z c (0,a) = f(a)-g'(-a) . (4.4) 

Combining ()4.3j) with (|4.4j) one finds {f'(x)} m = for m = j and {g'(x)} m = for m = k . 
Therefore the only possible final open strings are those of the form g(r — a) in the Z\ plane 
and fir + a) in the Z 2 plane. 

There is an exceptional case in which this is possible, namely when in the Z 1; Z 2 subspace 
there is a D brane that in the Zi, Z 2 subspace it appears as a D2 brane at an angle of ±7r/4 
with respect to the Z\, Z 2 planes. Using prime coordinates to denote the rotated system which 
is aligned with the D brane (located at, say, Z 2 < = 0), one has 

Z\i = cos 9 g{r — a) + sin 6 /(r + a) , Z 2 > = — sin 9 g{r — a) + cos 9 f(r + a) . (4.5) 

The coordinate Zy satisfies the M boundary conditions provided cos 9g{x) = sin 9f(x). The co- 
ordinate Z 2 i satisfies the T> boundary conditions provided cos 9f(x) = sin 9g(x) = sin 9 tan 9f(x), 
i.e. provided tan 2 9 = 1, implying 9 = ±7r/4. Note that in the rotated system the ring solution 
becomes 

Z v = V2L e iT cos a , Z 2 , = V2L e iT sin a . (4.6) 

Thus, with this D brane configuration at a special angle, the circular string can break into an 
open string. Successively, the open string can easily break in more fragments. This configuration 
is equivalent to a D brane such that on the Z±, Z 2 subspace it appears as a D4 brane with a 
flux Fi 3 = F 24 = 1. 

Quantum mechanically, a classically forbidden process can occur via tunnelling, but it will 
be exponentially suppressed for large mass (i.e. large L). To see this, we use the boundary 
state formalism. Omitting world-sheet fermions and ghosts, a general D brane is represented 
by the boundary state 

|B) = e -«^+4-4|o) , (4.7) 

where 0^,2? indicate the space coordinate mode operators with M or T> boundary conditions. 
We restrict our attention to the Z\, Z 2 subspace and consider a D-brane such that its projection 
on it is a D2 brane along the plane XyXy : 

\B) = e-^'^'-^'^'+^'^'+^'^l'lO) . (4.8) 
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Here all a operators refer to the world-sheet mode operators of frequency 1; the coordinate 
indices 1', 2', 3', 4' appear as subindices for clarity in the notation. Write 



ay = cos(#)ai + sm(9)a 3 , a 3 > = — sin(6 l )a 1 + cos(9)a 3 , (4.9) 
«2' = cos(#)a2 — sin(6 l )a4 , = sin(6*)a2 + cos(0)a4 . 

Then with our notation b^ ± = Ql ^ 2 , c± = a3 ^° 4 we get, 

_ e - sm(2e)(b\c\+J + b\+tf_^_+c^_)+^ ^ 

Therefore, the amplitude for the process in which the ring is "absorbed" by the brane thus 
becoming open strings on it, is 

A=(0\^±^\B) = (-sin(26)f . (4.11) 

Thus this process is exponentially suppressed unless 9 = ±7r/4. 

The rotating ring could still decay into open strings in a two-step process in which the ring 
first decays into a left-right symmetric (or antisymmetric) closed string state by emitting a 
closed string mode. Then the resulting (anti) symmetric closed string state can decay into two 
open strings. 

To calculate this process, we need to specify the emitted closed string state. Emission of 
massive particles from the ring are exponentially suppressed. One could still wonder if the 
process could take place by emission of a massless particle. We now show that this process is 
even more suppressed, like N~ 2N . 

We first consider an example, where we assume that the left-right symmetric closed string 
state is 

Then we need to evaluate the matrix element 

A R = ^ — =(0|(6c)^ : i.dX{l)e ik - XW : tf N \0) . (4.12) 

(#72)! VM lv ; v ; 17 v ; 

The calculation is straightforward and the final result is 

N ' /2 ( 9 -V 

A R = ^WM(-V2k)"-*-. (i /2S)* •£ r , (jv _ f ( + r T ) 1) , (f _ r) , + • ■ ■ («3) 

where dots stand for a similar term which is of the same order. To see the large N behavior, one 
can consider the two cases, N — N' = A r =fixed, or Affixed. Using the Stirling approximation, 
in both cases we find An = 0(N~ N ). One reason why this is so small is the factorial in the 
denominator (N — ^ + r — 1)! that starts from (y + ^ — 1)!, which is large. The only way 
that this would not be large is if the final state is similar to the ring state insofar as it contains 
b Nl and c Ni with both N — N\ and N — N2 small. But such states are far from being left-right 
symmetric or antisymmetric. They are small deformations of the rotating ring. 

One could consider more general emitted closed string states of the form (b'b') (c*&) 2 \0), 
N\ + N2 = N'. Since in this case we cannot have both N — N\ and N — N2 small, we again 
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expect a similar suppression. In particular, consider the extreme case (bW) N , i.e. N 2 = 0. 
This represents a decay of the rotating ring into the closed string state of maximum angular 
momentum, which, being a folded left-right symmetric string may successively decay into open 
strings by breaking. By a similar computation we find that the matrix element is of order 
0(N~ 2N ) for N — N' = JV =fixed and of order 0(N~ N ) for iV'=fixed. As expected for the 
reasons explained above, the process is suppressed. 

Thus we conclude that, except for a brane configuration at a special angle, the decay rate of 
the rotating ring in presence of D branes is still dominated by a decay into a similar ring state 
of lower mass by massless emission. The lifetime of the ring is therefore given by the formula 
of section 3.1, % ing = cog~ 2 M 5 . 



5 Case of compactified dimensions 

The results can be easily extended to the case in which some d c dimensions have been com- 
pactified. When some dimensions are compactified, type II superstring theories can have many 
exactly stable perturbative and non-perturbative states, associated with supersymmetric branes 
or supersymmetric fundamental strings wrapping cycles (in ten uncompact dimensions, only 
the DO branes of type IIA have finite mass). Here we are interested in the stability of non- 
supersymmetric configurations. We will focus on the rotating ring, since this seems to be the 
most stable non-supersymmetric (zero charge) state. We have explicitely considered the torus 
compactification, but the pattern is generic. We assume that the radii Ri of the compact dimen- 
sions are much smaller than the length L = \J~^- of our circular string (its mass is M = 

N being a large integer). In the opposite case, the circular string does not feel the effect of the 
compactification and we have the same results as in the uncompactified case. 

When two or more dimensions are compact, we have two different scenarios according to 
whether the compact dimensions occur in both the Z\ and Z 2 planes, where the circular string 
lies, or not. In fact, in the former case, say that X\ in the Z\ plane and X% in the Z 2 plane are 
compact, the circular string will wind around Xi with its Right part and around X 3 with its 
Left part, and it will easily break into winding modes. 

Instead, if the compact dimensions occur only in one of the above planes, say in Z\, then it is 
impossible for the circular string to break into winding modes in the Z\ plane while respecting 
the Virasoro constraints, since classically the Left part in the Z 2 plane would remain unbroken. 
Quantum mechanically, this process is expected to be suppressed exponentially for large length 
of the circular string. 

This has been numerically verified by extending the computation of [5] to the case of two 
compact dimensions [13]. 

Therefore, if the compact dimensions occur only in one of the planes, say Z\, the dominant 
decay channel will be by emission of radiation. 

It is easy to extend our operatorial computation of Section 3 to the case in which some d c 
dimensions are compactified. If no compact dimension occurs in the Z 12 planes, then the only 
modifications are the integration over the angles, because now we have to integrate J dQg-d c in 
the place of f dflg, and in the phase-space factor, which is now uj 

7-d cM ~2 R -l j n the pkce 

of lv 7 M~ 2 , because the momentum of the emitted massless particle becomes 9 — d c dimensional. 
In fact we can neglect the possibility of emitting a particle with Kaluza- Klein momentum as it 
would be effectively massive and thus exponentially suppressed. 
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The result is that the radiation spectrum is still described by a scaling expression 



2 S V 4 M 2 ^ I I I R^ ) f dc (a'Mu) . (5.1) 




The form of fd c (x) depends on d c , the number of compactified dimensions. The total radiation 
rate is obtained by summing over the integer No = |Mw getting 

r to t ~ g ^~ 6)/2 M-^ f[ . (5.2) 

i 

If some of the d c dimensions occur say in the Z\ plane, the pattern is the same as before, 
but the scaling function f$ c (x) is further modified by the fact that some among the momentum 
components k\ and &2, which appear in the expression giving fa c {x), are set to zero. If both 
the dimensions of the Z\ plane are compact and thus k\ — k^ = 0, we see from the expression 
of the Right contribution to the spectrum (jXTUj) , (|3.11|) that only N = 1 survives, that is, the 
spectrum is in this case a single point: 

uj = — — . (5.3) 
a'M v ; 

and the final state is similar to the initial one with the change M 2 = —N — ► Mf = —(N — 1). 

Consider, in particular, d c — 6 (corresponding to 3 + 1 uncompact dimensions). In the case 
the circular string lies on the uncompact dimensions in its, say, Right part and on a compact 
plane in its Left part, and its length is much larger than every compactification radius, it decays 

by emitting a massless particle of energy (J5.3j) and changing its length from <J to yj a 

This emission process takes a time T ~ -4tj Y\ 6 



The total time it takes to reduce its size from N large to N' of the order of some units, and 
thus easily breaking and disappearing, is obtained by integrating ^ ~ 1/T from N to 0. We 



obtain 



T „^ T . n _. (5 . 4) 



In a generic case in which some 6 — d sc radii are much larger than the circular string length 
L and d sc radii are of order 0(a'), the total time it takes for the string to decay completely in 
massless modes is 

( N o) 9s i V«' 

since now N will decrease on average by a finite amount (N ) of the order of unity in a single 
radiation step. 



6 Lifetime of average string state 

The lifetime of the average string state can be obtained by computing the decay rate averaged 
over all initial states of the same mass M. Here we will compute the average decay rate for 
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massless emission, and neglect the decay channel into two massive particles. Thus this will give 
an upper bound for the lifetime of the average string state. 

The average decay rate for massless emission was computed in [9] for the bosonic string. 
Here we generalize this calculation to the type II superstring. 

We are going to compute an amplitude of the form 

f , - M' 2 N 

? = const. g—a R a L , u> = = ^= , (6.1) 



where 



a R = -^i d 4z" N ' I ^-^[ii^yft*, l)*±flj*V R {k, 1)**] . (6.2) 



Nn r Ic z' Jc> z 2 

Vr is the vertex operator ()2.7|) and Mn r is the number of states at level N = N R = Nl in the 
NS-NS sector. Asymptotically, 

(O + l) 

Nn r = const. N R 4 e aVNR , a = ttVD - 2 , 
with D = 10 space-time dimensions. We introduce a basis of coherent states: 

oo oo 

|{An}> = n ex P[ A «-«n]l°;P> , IR}> = f[ exp [9 r ■ 1 0; p) , (6.3) 

n=l r=l— v 

where 9 r are Grassmann variables. Defining 



j a,a,6»,i 



lZ '. z '-N' I d± z -N ( g_ 4) 



•A/jv 7c 7c? 
x <{A„}, Wl-^^— Vj(A, 1)— ^—^(fc, l)z^|{A n }, {9 r }) , 

then the trace can be computed by integrating over all A n , A n and 9 r , 9 r with the weight factor 
e -9 r -e r e -\ n - A n ^ ue ^ norma li za tion. 

A shortcut to the calculation is to note that the following terms in (|6.2j) do not contribute 
to the trace: 

Y^{M^4^i4 k4\®N) = o , ^2{Q N \£4 k4\$N) = o • (6.5) 

The reason why these terms vanish is that they must be proportional to £.k or k 2 by Lorentz 
invariance. It is interesting to note that these correlators are not zero for particular states such 
as the rotating ring, but the whole expression vanishes after taking the trace. As a result, only 
the dX part of the vertex contributes. Then, since there is no fermion insertion, the fermion 
part in ()6.2|) gives the standard fermion partition function. What remains is the bosonic part 
computed in [9] plus the ghost contributions, which are implicit in ()6.2|) . 

In this case it is important to consider explicitely the ghosts, ensuring that we only keep 
physical states. In the be ghost sector, we compute the amplitude (Right part) 

*r= Yl (n N \c(i)\n NI )(n N ,\c(i)\n N } , (6.6) 

f2 jv j^jy' 
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where c(l) is the ghost field associated with the vertex operator V(k, 1) in (|2.7|) and |fijy) are 
"physical" ghost states (that is, such that the matter times ghost states satisfy the physical 
conditions). We then write: 

1 £\n Nf )(n N .\= <f ^z*- N, b C£mm)= <f ^~ N 'b t , (6.7) 

n N , J z * J z 

where the |\&) are all states built with b, c oscillators. The projector on states of level N' is 
the same of (|2.Hj) (here we use its ghost part), the bo insertion is due to the fact that while the 
sum is over all of the ghost states, the physical ones must satisfy the condition &o|^) = 0. 
These states, thanks to the b,c algebra, can be written as bo\^) where is any ghost state. 
Thus we have to compute 

A = Tr[(-1) F ° c{l)b c(z)z'*z*\ , (6.8) 

where F g anticommutes with all ghost modes. 

This can be done with the same coherent state formalism described above. The final result 

is 

00 

A = \\(l-w n )\ w = z'z. (6.9) 

n=l 

The superghost part of the average state computation gives directly the superghost partition 
function (since the vertex operators V(k, 1), V(k, z) do not carry any superghost contribution). 

The result is thus 2 

aR= LMpni^^i^.( 1 w -_ 1 (ir+8 , (I H , too) 



M N J z l-z N ~ N ' \y/z 
where 

00 00 00 00 

f(z) = l[(l-z n ) g 3 (z) = + z r ) , g A (z) = ]l(l-z r ) 92(z)=H(l + z r ). (6.11) 

n=l , — i , — i r=0 

For large N' the integral (|6.1Uj) can be computed by a saddle point approximation, the main 
contribution coming from z ~ 1. The computation is similar to the calculation of the asymptotic 
level density A/jv for level N. The behavior of the different functions in (|6.11|) at z ~ 1 follows 
from the modular properties of the theta functions. We find 

2 2 

f(z) = const. (1 - z) _ 5 e -B{f=ST j g 3 ( z ) ^ cons t. e TW^J ; (6.12) 



g A {z) = const, e ^5 ; g 2 ( z j ^ const, e 15 ^) . (6.13) 

The saddle point is at In z = with a = ~ 2 - Using N — N' = 2uy/N, a' = 4, we 

obtain for large N': 

ay/ N~2uVN-a\/N 

^M=C-^M 1 _ e _ au ^-^M — . (6.14) 

2 By the Jacobi identity, one can further simplify this expression setting ;4y<7§ — T/jflf = <?i- 
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Thus, by eq. (|6.1|) . the average decay amplitude is given by 

f = const. g 2 s £ 2 i 2 cr grc (w,T) -J— oj 8 , (6.15) 

cr — 1 

where T = is the Hagedorn temperature of type II superstring theory and a gvcy {uj) is a 
grey body factor given by 

1 - — 
1 — e t 

a grey (^,T) = uj- . (6.16) 

(1 — e 2T ){1 — e 2Tj 

Thus an average superstring state in type II superstring theory emits radiation as a grey body 
at the Hagedorn temperature. In the case of type I open superstrings, the spectrum for massless 
vector emission is obtained from the right part (J6.14)) . i.e. 

f = £ • £ M~ l 6 T \ u 8 , (6.17) 
1 - e T i 

where T\ = is the Hagedorn temperature of type I superstring theory. Hence average 
states in type I superstrings emit radiation as a perfect black body. 

To compute the total average decay rate, one sums over all Nq. Recalling that uj = Nq/2M, 
this sum can be approximated by an integral over u by introducing the differential 2Mduj. We 
find 

f tot ~^ 2 M. (6.18) 

Since T represents the average decay rate taking into account massless emission only, the lifetime 
of an average state has the upper bound 

T<gfM- 1 . (6.19) 

The same result holds for the bosonic and heterotic string theory. 3 



7 Discussion 

The average lifetime computed here, T < M _1 , shows that most of the massive states at each 
mass level have a short lifetime, which decreases for larger masses. There can only be a few 
states at each mass level whose lifetime increases with the mass. Some of these few states are 
given by the family of states \^ n ,k)i N = 2k, studied in [5], interpolating between the rotating 
ring n = k and the state n = 0. The corresponding classical configurations describe rotating 
elliptical strings, which become circular for n = k and a folded, straight closed string for n = 0. 
For generic n one finds that the lifetime T also increases with the mass, though at much lower 
rate for states that significantly differ from the rotating ring (i.e. states with n not close to k). 
When n k, the ellipse is very eccentrical and the decay by breaking into two massive string 
states becomes the dominant channel. These are all states of high angular momentum. 

The law T ring = const. M 5 has been derived in four independent ways. 1) By evaluation of 
the imaginary part of the two-point function [5] , 2) by operator calculation of the full massless 

3 Sincc in open string theory the decay in massive states is important, because breaking is never suppressed, 
the bound put by the massless emission result 7ty pc i < const. \fa! in this case could be rather weak. 
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emission rate, 3) by operator calculation of the dominant decay channel into a ring of lower 
mass and 4) by classical radiation from the rotating string soliton configuration. The latter 
also confirms the identification of the quantum state with the string soliton. 

The rotating ring has the special feature that it decays predominantly into another rotat- 
ing ring, i.e. the decay can be viewed as a rotating circular string whose radius gradually 
decreases as it emits radiation. The decay rate of the rotating ring has the scaling property 
T = eg 2 M~( 5_dc ) f(N ), where N = M 2 — M' 2 = 2uM, d c = number of compact dimensions 
with radius R « a'M, and a spectrum sharing some qualitative features with a thermal spec- 
trum, such as an exponential tail. This is not the case e.g. for the state of maximum angular 
momentum (tftf) N \0). 

In the case of the average state, one finds as in [9] that the spectrum becomes exactly 
thermal in the large mass limit and that the average type II superstring emits radiation as 
a grey body with temperature equal to the Hagedorn temperature (while the average type I 
superstring behaves as a perfect black body at the Hagedorn temperature). 
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Appendix A Notation 

We consider type II superstring theory in the NS formulation. The world-sheet action is given 



by 




(A.l) 



We will be interested in string states corresponding to excitations in two planes, 



X 1 + iX 2 2 _ X 3 + iX 4 
V2 ' y/2 



v/2 y/2 



The expansion for the bosonic world sheet fields are as follows 
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\Pni: i b m± ] 3nm j [^n± j ^m± ] ^nm j [Cn± > ^mi ] S nm , [c n -|- , C m _j_ ] 5 n 

For the fermion fields, we have 



r r 

The mass formula is given by 

a'M 2 = 2(N R + N L ) , N R = N L , (A.3) 
where the number operators are 

oo oo 
Nr = ^ n ( b l+ b n+ + bl_b n _ + c{ + C n+ + c{_C n _) + ^ *i>-W? ~ a • ( AA ) 



n=l i — I 

' — 2 



The expression for is similar, with the obvious changes. The normal ordering constant a 
is a = 1/2 in the NS sector and a$ = in the R sector. 
The basic correlators are 

(d w X(w)d z X(z)) = -^j-—r 2 , mwMz)) = ^773V • ( A - 5 ) 
We use the following angular coordinates 

x\ = sin 9 cos a cos /? cos <p , x 2 = sin 9 cos a cos /3 sin ip , 

x 3 = sin 9 cos a sin (3 cos -0 , £4 = sin 9 cos a sin (3 sin z/> , 
£5 = sin # sin a cos /3' cos <// , xq = sin sin a cos f3' sin <// , 
£7 = sin 6 1 sin a sin cos ■?/>', xs = sin sin a sin sin ?// , x$ = cos # . 
Thus the volume element is 

<i 8 f2 = sin 7 9 cos 3 a sin 3 a cos j3 sin /3 cos /?' sin (3' d9dadpdtpdipdp' 'dip' dip' 

In order to sum over polarizations, we use the relations 

hi Uj 

= E^ = rf - 1 = 8 - ( A - 6 ) 

This gives (k and £ are defined in (|3.5|L ()3.6|) ) 

^ IZJ2 ^ IJ.I2 ^ IL.I2 v ' 

P I I £ I I £ I I 

Writing 

AT 2 

o 2 iv 2iip 2 2 _ • 2 /) 2 2 a 

2k = — — e v v , v = sin cos a cos a , 
4iV 



we have 



£ £ £ 
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Appendix B Other contributions to the decay rate 

of the rotating ring 

The remaining contributions to the decay rate of the rotating ring indicated in section 3.1 are 
a 2 = -i-L l d A ^(0|6 w e- ifc - x(1) e ifc - xw e9X(l)6 tW |0)(0|^^W^^(^)^!i|0) , 

N\ J Z 2 2 

a 3 = i_L i^L ^( |^ e - ifc - x(1) e ifc - x( ^.9X(z)6 tW |0)(0|V'fe^(l)A:^(l)^!i|0) , 

N\ J Z 2 2 

a 4 = --L /— z^(0|6^e- ifc - x(1) e ifc - x(2) 6 tW |0)(0|^fe^(l)^^(l) ^Mz)kMz)ip\ |0) , 

N\ J Z 2 2 

a5 = J_ ^(0|6 Jv e- ifc - x «e 4fc - x ^6 tJV |0)(0|e9X(l)e.9X( 2 )|0)(0|^ z 1 |0) 

N\ J Z 2 2 

Computing the matrix elements, we find 

t fdz N n o o N i , iV!(2KK) m (l-z) 2m+1 

^ m=0 v 7 



0" 5 



J ^ m\ 2 (N - m)\ z m 



l-z)_ 



and (J3 = <j\. Expanding the binomials (1 — z) and computing the integral over z, we find 

'-l) m ~ Na N\{2m + l)\(2KR) m [(eeRK-e 2 R 2 )(m -N +l)-(e 2 K 2 -eeRK)(m + N + 1 



2 



m!(m + 1)!(JV - m - l)!(m - iV + l)!(m + iV"o + 1)! 



^ (-l)^ojv!(2m)!(2««)» 

a 4 = 4£ • £ KK 2^ 



r m! 2 (iV-m)!(m-iV"o)!(m + iV"o)! 



ot t ^ (-l) m - JV °iV!(2m-2)! (2««) 



: m! 2 (iV - m)!(m + iV - l)!(m - JV - 1)! ' 

Now we can estimate which is the dominant contribution in the large N limit, with fixed N — N'. 
Using Stirling formula (J3.24J1 and the fact that kk ~ iV -1 , we see that the dominant term is a±, 
with 0i = O(N). Similary, we see that 02,03,05 are of order 0(1) and 04 is of order 0(iV _1 ) 
(the explicit numerical evaluation shows that 02, •••,05 are two or more orders of magnitudes 
smaller than <j\ already for N = 20). The application of the Stirling formula is justified in the 
present case because decays with N ^> 1 are always suppressed as long as N is large. 



22 



Appendix C Tail of the spectrum 

We consider the expression for En f)3.18|) . Using the Stirling formula (|3.24|) . we find 



_ x y+s 



E u = ^N 2 x 2N ° y — ^a(r)a(s) , (C.l) 

V ^ T(2N + r + s + ^ V 1 K ' 



a(r) 



r,s=0 11 2 > 

(2r + 2N )\ 
~~ r\(N - 1 + r)\(r + 2N ) 



TV 

where x = 2kkN = We are interested in the large N behavior of this sum. Because the 
series converges very rapidly, we can approximate 

(2N + r)\ - {2No) ' (2iV + S )! " {2No) " (C ' 2) 
At the same time, we use Stirling formula to approximate 

r(2iV + r + S + y) - {2N )\ (2N ) r+s+ % 

V2^(2N ) 2No+r+s+6 e~ 2N ° . (C.3) 



We find 



2N 2N 

- a "^wsr** s{Nof • (a4) 



We can write this sum as 



S(N ) = ^~jj t N °-' e*p[-xt + x = ^ . (C.6) 

The integral can be computed by saddle-point approximation, with the result 

S(N ) = const. N$ e ^iog2-JV iogjV ? ( C ^ 

where (3 is a constant of order one. Adding the phase space factor, we finally obtain 

T n = -i^/£n = N( J^) 7 Arf-6 e - 2 7Vo(2io g2 -i) (C-8) 
M 2 K 2^N J 

i.e. 

r u = TV"! Nl +2P e -2JVo(21og2-l) (C 9) 

Similarly, one can consider E22 and E12. We find the same behavior. 
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